Spin-Glass and Chiral-Glass Transitions in a ± J Heisenberg Spin-Glass Model in 

Three Dimensions 
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The three-dimensional ±J Heisenberg spin-glass model is investigated by the non-equilibrium 
relaxation method from the paramagnetic state. Finite-size effects in the non-equilibrium relax- 
ation are analyzed, and the relaxation functions of the spin-glass susceptibility and the chiral-glass 
susceptibility in the infinite-size system are obtained. The finite-time scaling analysis gives the 
spin-glass transition at T Bg / J = 0.21^q'q2 an d the chiral-glass transition at T cg /J = 0.22^Q'gg. The 
results suggest that both transitions occur simultaneously. The critical exponent of the spin-glass 
susceptibility is estimated as 7 sg = 1.7 ± 0.3, which makes an agreement with the experiments of 
the insulating and the canonical spin-glass materials. 
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It has been an unsettled issue whether a short-ranged 
spin-glass model can explain the spin-glass transition ob- 
served in real materials. |l| The existence of the spin-glass 
transition in the Ising model has been widely recognized: 
various types of investigations [^|j3| gave consistent values 
of the transition temperature. The recent large-scale sim- 
ulation study Q evaluated the critical exponents which 
quantitatively agree with the measurements of the Ising- 
type spin-glass material, Feo.5Mn . 5 Ti03. pj However, 
the spins and the magnetic interactions of many spin- 
glass materials are mostly isotropic, and the anisotropy 
effect can be regarded as a secondary effect. Therefore, 
the Heisenberg model with random interactions is consid- 
ered to be the simplest theoretical model. The problem 
is whether this model exhibits the spin-glass transition or 
not. Previous numerical investigations suggest that the 
finite-temperature transition does not occur. |5|-|l 1[| 



In order to solve this contradiction, Kawamura |7j] in- 
troduced the chirality mechanism. The chirality is a 
variable whose absolute value corresponds to a volume 
of a parallelepiped spanned by three spins, and whose 
sign represents the handedness of the spin structure. In 
the chirality mechanism, the chirality is considered to 
freeze randomly at a finite temperature without the con- 
ventional spin- glass order. The experimentally-observed 
spin-glass anomaly is considered as an outcome of a mix- 
ture of the chirality with the spin by the weak magnetic 
anisotropics. The existence of the chiral-glass transi- 
tion has been confirmed by several numerical analyses. 
Jf|— |i~0[| However, it was pointed out recently by one of the 
present authors (S. E.) that the finite-temperature spin- 
glass transition may exist by examining the stiffness of 
the system when all the spins on a surface are rotated and 
reversed, respectively. fl^,|l3| The same evidence was also 
reported by a dynamical Monte Carlo simulation which 
removes the global spin rotation. 

In this Letter, we focus on the existence or the non- 
existence of the spin-glass and the chiral-glass transitions. 



By employing the non-equilibrium relaxation method, we 
have revealed that both transitions may occur simulta- 
neously at a finite temperature T/J ~ 0.21. This re- 
sult possibly settles the discrepancies between the exper- 
iments and the theories in the Heisenberg spin glass in 
three dimensions. 

In the simulational studies of the phase transitions, one 
obtains equilibrium quantities in finite-size systems, and 
the finite-size scaling is utilized to extract the thermody- 
namic properties from them. However, the dynamics of 
the simulations become very slow due to the randomness 
and frustration in the spin-glass models, which makes 
it hard to realize the equilibrium state in large systems. 
The system size accessible by the equilibrium simulations 
is restricted to linear sizes L ~ 30 up to now. In this sit- 
uation, reliability of the finite-size scaling might become 
doubtful, which causes controversial conclusions of the 
existence [[t5| or the nonexistence raj of the spin-glass 
transition depending on the sizes and the temperatures 
of the data used. In this Letter, we follow a completely 
alternative approach to the thermodynamic limit, i.e., 
we observe the relaxation of the infinite-size system to 
the equilibrium state. In order to extract the equilib- 
rium properties, the finite-time scaling analysis is utilized 
instead of the conventional finite-size scaling. This ap- 
proach is known as the non-equilibrium relaxation (NER) 
method. Jl^-Q Actually, we prepare a very large lat- 
tice with the paramagnetic state ( perfectly random spin 
configuration) , and observe the relaxation of the suscep- 
tibility. The simulations are stopped before the finite-size 
effect appears. Time dependences of the susceptibility di- 
verge algebraically as %(t) ~ t^^ zv JT^-pl]] in the critical 
region, while it converges to a finite value in the para- 
magnetic region. We detect the phase transition by this 
difference. Another advantage of this strategy is an ab- 
sence of a trivial exponential decay due to the global spin 
rotation. Since the spin autocorrelation function is essen- 
tially a one-body function about the spin, one may en- 
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counter this exponential decay even in the ordered phase 
after a time scale of the global spin rotation. E3] 

We consider the ± J Heisenberg spin-glass model on a 
simple cubic lattice of TV = Lx Lx (L + l) under the skew 
periodic boundary conditions, Ji = Jij^i ' Sj- 

The exchange interactions Jij take a value of + J or — J 
with the same probability, and Si is a three-component 
vector spin with \Si\ = 1. The summation runs over all 
the nearest-neighbor spin pairs. We evaluate the spin- 
glass susceptibility, x sg , and the chiral-glass susceptibil- 
ity, x C g, through the following relations. 

^ = ^K(E c S0 2 >]- ( 2 ) 

where q^. v — (1/iV) J2i S^Sf'l is a replica overlap be- 
tween the /i component of a spin i on a replica (a), S^ }, 

and the v component of the spin on a replica (6), sfl- 
The bracket [• • -] c denotes the configurational average, 
while (• • •) denotes the thermal average. The indices, // 
and v, stand for three components of spins: x, y, and z. 
The local chirality C z -° is defined by three neighboring 
spins as (M = S (a l ■ {S\ a) x S {a) - ), and e M denotes 

a unit lattice vector along the /i axis. 

We use the single-spin-flip algorithm with the heat- 
bath probability, fl For a given random bond config- 
uration, eight replicas are prepared with different ran- 
dom initial spin configurations. Each replica is updated 
in parallel by a different random number sequence, and 
twenty eight replica-overlaps are calculated and averaged. 
Note that these data are correlated with each other, and 
only seven among them are independent. Then, we take 
averages of Xsg(i) and Xcg(t) over the bond distributions. 
Typical numbers of the bond configurations are tens to 
thousands depending on the sizes, the temperatures, and 
the time ranges. 

From the non-equilibrium relaxation functions of x sg 
and Xcg, we extract the equilibrium properties and esti- 
mate the transition temperature by the finite-time scal- 
ing analysis, |2^,^4| which is a direct interpretation of 
the conventional finite-size scaling through the relation, 
r ~ £ z . We determine the correlation time r(T) at 
each temperature and a combined exponent A(= j/zu) so 
that the scaled functions, for various tempera- 

tures fall onto a single curve when plotted against t/r(T). 
Then, we extrapolate t(T) with an assumption that they 
diverge algebraically when the temperature approaches 
the transition temperature T c as r(T) oc (T — T c )~ zv . 
We obtain the most probable estimate for T c by the least- 
square fitting of this form. 

In the finite-time scaling, it is crucial to use the data 
free from the finite-size effect. We estimate a characteris- 



tic time that the finite-size effect appears by changing the 
system size from L = 9 to L = 59 at each temperature. 
The typical example is shown in Fig. [I]. The temper- 
ature is within our estimate of the spin-glass transition 
temperature obtained afterwards. The finite-size effect 
of L = 19 data already appears at t = 1000 steps, af- 
ter which a long relaxation to the finite-size equilibrium 
state is waiting. This "finite-size" characteristic time is 
roughly scaled as r ~ L z with the dynamic exponent 
z ~ 5. By this relation, we confirm that the relaxation of 
L — 59 before t = 100000 steps is not influenced by the 
finite-size effect, and thus can be considered as that of 
the infinite system. At each temperature, we check this 
time scale and use only the data with the approval. 
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FIG. 1. The size dependences of x S g(i) (up to L = 59) and 
Xcg(i) multiplied by 100 (up to L = 29) at the temperature 
T/J — 0.21. After the initial relaxation, XsgW an d Xcg(£) 
converge to critical behaviors t x (dotted straight lines) with 
A = 0.35(sg), and A = 0.23(cg). These exponents coincide 
with the estimates obtained by the finite-time scaling. 

In this figure, we notice that the appearance of the 
finite-size effect is in opposite ways between the spin-glass 
susceptibility and the chiral-glass susceptibility. The for- 
mer is always underestimated, while the latter is over- 
estimated showing upward bending. This tendency con- 
tinues to the equilibrium values; the chiral-glass suscep- 
tibility always decreases as the system size is enlarged. 
This size dependence is quite strange as the susceptibil- 
ity, but it can be understood by considering a definition 
of the chirality. Since the chirality corresponds to a vol- 
ume of a parallelepiped spanned by three neighboring 
spins, it becomes larger when the relative angles between 
the neighboring spins are enhanced by the finiteness of 
the lattices. Thus, we must pay much attention to the 
finite-size effects in the relaxation of the chirality. For 
example, at paramagnetic temperatures, the chiral-glass 
susceptibility converges to a finite value in the large size 
limit, while it shows an upward bending in the small size 
limit. There exists a crossover size in which the non- 
equilibrium relaxation looks algebraically divergent as if 
it exhibits the criticality. This crossover size is consid- 
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ered to correspond to the chiral correlation length. If 
one tries to detect T cg solely by finding the algebraic 
relaxation behavior, he may be misled to a wrong tem- 
perature where £ = L instead of £ = oo. Therefore, the 
chiral-glass transition temperature tends to be overesti- 
mated, and the spin-glass transition temperature is to be 
underestimated, if the system size is insufficient. 

As shown in Fig. [l], the initial relaxation of the spin- 
glass susceptibility continues until t ~ 10 3 steps, then it 
exhibits the critical relaxation. In the finite-time scal- 
ing analysis, we must discard this initial relaxation data. 
It is also noticed that the finite-size effects of the small 
sizes {L < 19) appear before the initial relaxation ends. 
The critical phenomena may not be extracted from the 
equilibrium states of such small systems. 
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FIG. 2. (a) The finite-time scaling of Xsg(i) with 
A = j/zi/ — 0.36 and Xcg(t) (multiplied by 300) with A = 0.24. 
The t _A are depicted by straight lines. The data of L = 59 
at T/J =0.32, 0.29, 0.28, 0.27, 0.26, 0.25, and 0.24 are plot- 
ted, (b) Extrapolations of the correlation time t(T) of Xsg 
and Xc g obtained by the scaling. The transition temperatures 
are estimated so as to minimize the least-square errors. Six 
temperatures ranging 0.24 < T/J < 0.29 are used. 

Results of the finite-time scaling of the spin-glass sus- 
ceptibility and the chiral-glass susceptibility are shown 
in Fig. ^. The combined exponent A and the correlation 
time t(T) are estimated by using the data of L = 59 
until t = 70000 steps at seven temperatures ranging 



0.24 < T/J < 0.32. The data in the initial relaxation 
are discarded. The scaling function converges to a t~ x 
line in the equilibrium region (t/r(T) — > oo), where the 
susceptibility converges to an equilibrium value in the in- 
finite system. On the other hand, it becomes a flat line in 
the non-equilibrium limit, t/r(T) — > 0, since the suscep- 
tibility behaves as %(i) ~ £ A in the critical region. There- 
fore, this curve interpolates the non-equilibrium and the 
equilibrium limits. 

The correlation time is extrapolated by (T — T^)~ zv , 
and we obtain the transition temperature that gives the 
best least-square fitting. The exponent is estimated to be 
zv = 4.7 ± 0.8. Combining this value with A = 7 / ' zv — 
0.36 ±0.02, the exponent for the spin-glass susceptibility 
is estimated. 



T / 7 — 21 +001 



7sg = 1.7 ±0.3. 



This transition temperature agrees well with that ob- 
tained by the defect free-energy method, fl3| l and that 
obtained by the dynamical simulations. |L4j The expo- 
nent is consistent to that of the non-linear susceptibility 
in real spin-glass compounds: 7 = 2.2 ± 0.1 (CuMn & 
AgMn) j§ and 7 = 2.3 ± 0.4 (Cd Cr2xo.s5In2xo.85S4). 
p6"t This may suggest that the present theoretical model 
is relevant to the spin-glass materials. Estimates and the 
comparison of other exponents are left for a future inves- 
tigation. It requires more detailed and accurate calcula- 
tions of the Binder ratio, which diverges in the critical 
region as t d / z . Preliminary calculations of the Binder 
ratio, however, suggest that the size effects and the sam- 
ple dependences are much severer than those of the sus- 
ceptibility. The algebraic divergence is only observed in 
lattices larger than L = 19, below which the finite-size 
effects appear in the initial relaxation. The details will 
be reported elsewhere. 

Analyses of the chiral-glass transition are done in the 
same procedure. Note that the amplitude of % cg is 
smaller than that of % sg about two digits, which causes 
rather large numerical errors. A possible value of A 
ranges as A = 0.23 ± 0.04. The transition temperature 
and its exponent are obtained as 



T cg / J = 0.22; 



-0.01 

-0.04' 



1.0 ±0.5. 



Note that this value is consistent with the estimate ob- 
tained by the equilibrium simulations which gives T cg ~ 
0.21. [B7| The value of 7 cg does not contradict with that 
of the model with the Gaussian bond distributions. JTcj ] 
The exponents A for \ sg and \cg estimated in the scaling 
also coincide with the slope of the relaxation functions 
at the transition temperature in Fig. ^. This evidence 
supports the validity of the present analysis. By the re- 
sults presented in this Letter, it is conjectured that the 
spin-glass transition and the chiral-glass transition occur 
at the same temperature. Since the chirality must freeze 
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if the spin freezes, the transition is possibly driven by the 
spin degrees of freedom. 

We have observed the simultaneous finite-temperature 
spin-glass transition and the chiral-glass transition in the 
±J Heisenberg model in three dimensions. This was 
made possible by observing the non-equilibrium relax- 
ation of the infinite-size system started from the para- 
magnetic state. We consider the inconsistency of the 
present results with those of the previous investigations 
[p|-|ri| can be explained by the differences of the lattice 
size and the temperature range investigated. Generally, 
one expects to observe the critical behavior as long as the 
lattice size is comparable with or larger than the correla- 
tion length, L > £. Contrary, when the correlation length 
exceeds the size, L <C £, a very large but a finite cluster 
swallows the whole system, and only a trivial size effect is 
observed. As the temperature is lowered from the para- 
magnetic phase, the correlation length exceeds a lattice 
size at a temperature above the transition temperature. 
If one mixes the data above and below this crossover 
temperature, the finite-size scaling analysis may mislead 
to a wrong conclusion. Generally, frustration causes a 
longer correlation length, and the size of the system is 
sometimes essential to detect the phase transition par- 
ticularly in frustrated systems. p4| In the present non- 
equilibrium relaxation analysis, we have used only data 
of the temperatures and the time ranges where they can 
be considered as the infinite system. The time range is 
limited within 10 5 steps, however, it is enough to equi- 
librate the system at higher temperatures, T/J ~ 0.3, 
where we start the scaling analysis. We gradually lower 
the temperature so that the scaled relaxation functions 
overlap with each other to ride on a single curve. At 
lower temperatures, T/J ~ 0.24, 10 5 steps correspond to 
a time scale that the finiteness of the correlation length 
begins to appear, i.e., the relaxation function begins to 
bend. In order to see a longer relaxation, we need to 
enlarge the system size to eliminate the size effect. This 
must become much more time-consuming and we leave it 
for the future study. However, it should be noticed that 
the finite-size effect appears quite early in the relaxation 
process as shown in Fig. [lj and what is slow is the re- 
laxation to the equilibrium state after this "finite-size" 
characteristic time scale. The relaxation of the infinite 
system might be simple: after an initial relaxation of 
10 3 steps, the critical relaxation continues to infinity. As 
for the relation between the model and real spin-glass 
materials, the critical exponent obtained by the scaling, 
7s g ~ 1.7, agrees with the experiments. 
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